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Abstract. Let B n denote the nth Bell number. We use well- 
known recursive expressions for B n to give a generalizing recursion 
that can be used to prove Touchard's congruence. 



1. Introduction 

For a positive integer n, the Bell number B n is the number of ways 
a set of n elements can be partitioned into nonempty subsets. Com- 
putations of Bell numbers rely on well known recursive formulae. For 
example, one can compute B n+ i by enumerating partitions accord- 
ing to the size of the subset S which contains the n + 1st element; if 
|«S| = n + 1 — d then there are ( choices for the other elements of 
S and Bd ways to partition the elements not in S. Hence we have 



'1-1) B n+ i = B a 



d=0 



n 

n — d 



Alternatively, one can compute B n by keeping track of the number of 
subsets in a given partition; the Stirling number of the second kind 
counts the number of partitions of n elements into precisely k subsets, 
and so one has the well-known formula 



;i-2) B n = 



k=l 



Calculations for {^} rely on the binomial-like recurrence relation 



Using the identities above, we derive an expression for B n+ j which 
generalizes both (11. ip and (11.21) : 
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Theorem 1.4. For positive integers n and j , 

k=i ^ * 

where Pj(x) is the degree j polynomial 5]J =0 5j_ r g3; r '. 

It appears that this formulation wasn't discovered until recently 
when a combinatorial proof of this result was given in [I]. We arrived 
at this result independently through algebraic manipulations of (11.11 - 
11.31) ; since this result comes directly from these well-known identities, 
it is surprising that Theorem II .41 wasn't recognized much sooner. 

This project began when the first author was a student in the second 
author's elementary number theory class. Before knowing the general 
form of the polynomials Pj(k) from Theorem 11.41 the authors could 
find specific identities such as 

n 

B n+5 = J2( k5 + 5fc4 + 20fc3 + 50fc2 + 7 $k + 52) 
k=i 

These computations were being carried out just as modular arithmetic 
was being introduced in the concurrent course, and it was clear that 
when j was prime the polynomials Pj(k) were ripe for simplification 
modulo j. Flushing this observation out not only provided a great tour 
of some of the most familiar identities and techniques for computations 
modulo p, but ultimately led to a rediscovery of 

Corollary 1.5 (Touchard's Congruence). For positive integers n and 
m and p a prime number, 

B n+pm = mB n + B n+ i mod p. 



2. A Proof of Theorem 11.41 

We prove the result by induction, with j = our (trivial) base case. 
By induction we have 

B„ +1+1 = £ Pj (k) { " + } = £ £ Bj _ P ) k r 

k=l ^ ' k=\ r=0 ^ ' 
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Applying identity (Oil to j™^ 1 }, and using {"} = {J^} = 0, we have 



n+l j 



k=l r=Q 



r I VIA; — 1 1 1 A; 



n+l j / \ ( ^ n j 



fc=2 r=0 x / ^ fc=l r=0 



r / 1 /c 



= E E b,- (j) <* + irH + £ t (jV« 

fe=l r=0 V ' ^ J k=l r=0 V ' 

The coefficient of i n this sum is then 

e^C)0 +b -"-"G-i) 



We now make the change of variable d = j — r in the first sum and 
apply identity (II. ip . leaving us with: 



vd=0 v ^ 



Thus we have B n+j+1 = ]T E B ^ ( J + ^ ^{l} = 2 ^i+iW It} 



k=l r=0 x / k j k=1 



3. Computations modulo p 

If p is an odd prime, it is well known that ( p r ) = mod p whenever 
< r < p m , and so all but two of the terms of P p ™(x) are congruent to 
zero modulo p. Applying Theorem 11.41 and Fermat's Little Theorem, 
we therefore have 

B n+pm = J2(B pm + k pm B ) j" j = J2(B pm + k) i n \ mod p. 

k=l ^ ' k=l ^ ' 



4 



G. HURST AND A. SCHULTZ 



Since P x (k) = l + h, Theorem O gives B n+1 = £(1 + A;){™}, and so 
rearranging the previous congruence gives 

B n+pm = (B pm - 1) \l \ + J> + *) 
k=i ^ ' k=\ 
= (B p m - l)B n + B n+1 mod p. 

The same congruence holds for p = 2 as well: if m > 2 then ( 2 r ) = 
mod 2 for < r < 2 m and our previous argument holds, and if m = 2 
the only additional term is 3/c 2 i?2 = Qk 2 = mod 2. 

To verify Corollary ll.5[ then, we only need to prove the following 

Lemma 3.1. For every positive integer m and prime number p, B p m = 
m + 1 mod p. 

Proof. B p m enumerates the partitions of Z/p m Z. Our strategy will be 
to let Z/p m Z act on these partitions in the natural way: for elements 
x and y of Z/p m Z we define f y (x) = x + y mod p m , and f y (P) is the 
partition we get by applying f y element-wise to P. Any partition not 
fixed under this action will belong to an orbit of size a power of p, and 
so the number of fixed partitions is equivalent to B p m modulo p. 

So suppose you have some fixed partition P with elements a and b 
inside subsets A and B (respectively). Then clearly f^ a (a) = b, and 
since P is fixed this means fb- a {-A) = B. Hence for a fixed partition 
P, all subsets of P must be the same size, and therefore some power 
of p. We claim that the only fixed partition whose subsets are size p 1 
is the partion whose subsets contain elements which are congruent to 
each other modulo p m ~K This would leave us with m + 1 many fixed 
partitions, as desired. 

To prove the claim, suppose to the contrary that we have a fixed 
partition P with elements a, b in the same p J -element subset A which 
satisfy a ^ b \Tiodp m ~K Now fb-a( a ) = b, and so A is permuted 
by the action of fb- a - Hence for any integer r, the r-fold composition 
of the map fb- a — namely, the map f r (b- a ) — again takes a to some 
element of A. Now clearly f r (b-a)(a) = fs(b-a)(a) mod p m if and only 
if 

a + r{b — a) = a + s(b — a) mod p m . 
Since a ^ b mod p m ~\ however, this congruence forces r — s = 
mod pi +1 . Hence for r between 1 and p* +1 , the elements f r (b-a){.o) are 
distinct contituents of A. It follows that |*4| > p> , a contradiction. □ 
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